Entanglement spectrum of the two dimensional Bose-Hubbard model 
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We study the entanglement spectrum (ES) of the Bose-Hubbard model on the two dimensional square lattice 
at unit filling, both in the Mott insulating and in the superfluid phase. In the Mott phase, we demonstrate that the 
ES is dominated by the physics at the boundary between the two subsystems. On top of the boundary-local (per- 
turbative) structure, the ES exhibits substructures arising from one-dimensional dispersions along the boundary. 
In the superfluid phase, the structure of the ES is qualitatively different, and reflects the spontaneously broken 
U{1) symmetry of the phase. We attribute the basic low-lying structure to a so-called "tower of states" (TOS) 
Hamiltonian of the model. We then discuss how these characteristic structures evolve across the superfluid to 
Mott insulator transition and their influence on the behavior of the entanglement entropies. Finally, we briefly 
outline the implications of the ES structure on the efficiency of matrix-product-state based algorithms in two 
dimensions. 



Introduction — In recent years the cross fertilization be- 
tween quantum information and condensed matter has led to 
new insights into the physics of low dimensional systems [1]. 
In particular the concept of entanglement spectrum has estab- 
lished itself as an informative and intriguing theme. Consid- 
ering a bipartion of a system into parts A and B, the entangle- 
ment spectrum (ES), {^i}, is defined in terms of the Schmidt 
decomposition 

|^)=^e-«'/2|V.^)0|Vf). (1) 

i 

Here |?/;) is the ground state, and the states j-i/)^) (IV'i^)) form 
an orthonormal basis for subsystem A (B). The ES ~ 
— logAi} can also be thought of as the spectrum of the so 
called entanglement hamiltonian He = — log pa where the 
reduced density matrix pA is obtained after tracing out the B 
part of the system density matrix The ES can also 

be used to construct entanglement entropies (Renyi and von 
Neumann) which quantify the entanglement between the two 
subsystems. 

Many results are now available for the ES and entangle- 
ment entropies of one-dimensional (ID) systems. In contrast, 
higher dimensions are far less explored. Much of the ES lit- 
erature on 2D systems focuses on topological phases [2]. A 
detailed understanding of generic ES features for more com- 
mon 2D systems is not cuiTently available. 

In this work we present a thorough investigation of the ES 
for the 2D Bose-Hubbard model [3, 4]. We extract features 
for both the gapped Mott insulator and the superfluid (gap- 
less, symmetry-broken) phases. This analysis provides in- 
sights into the ES which, in addition to the intrinsic impor- 
tance of the Bose-Hubbard model, are likely to be generic to 
gapped and symmetry-broken gapless 2D phases. Our analy- 
sis is based on DMRG calculations on a cylindrical geometry, 
complemented by perturbative calculations (Mott phase) and 
by analysis of the tower of states analogy (superfluid phase). 

In the Mott insulator, the ES is "boundary-local" in the 



sense introduced in Ref. [5] for ID gapped systems. In ID 
gapped phases adiabatically connected to a "simple" state 
(product state or low-entanglement state), the ES is organized 
in a perturbative or hierarchical structure closely linked to 
boundary-locality: higher levels in the ES coiTesponding to 
excitations farther from the boundary [5]. The present work 
shows how ES features of the 2D Mott insulator can be under- 
stood from this viewpoint. In 2D the boundary is an extended 
object and not a single point. In the ES, this shows up as 
a multiplicity structure superposed on top of the perturbative 
one: there are multiple ES states for each order of the per- 
turbative hierarchy. This reflects the physics that excitations 
created close to the boundary can move without changing the 
distance from the boundary, resulting in dispersion structures 
within the ES at every order This can be cast in the lan- 
guage of entanglement Hamiltonians ("H e) for which we pro- 
vide some analytic perturbative expressions. We also provide 
analytic perturbative expressions for several entanglement en- 
tropies (von Neumann and Renyi) in the gapped phase, and 
discuss how area laws emerge in the perturbative framework. 
We believe that these features of the ES in the Mott phase 
(boundary-locality, multiplet structure, ID dispersions, etc) 
are general for (non-topological) 2D gapped systems (gapped 
fermionic models, Heisenberg XXZ, etc.). 

At the transition between the Mott and superfluid phases the 
ES structure changes dramatically. We show that its lower part 
(which we call ES "envelope") sharply reflects the different 
nature of the two phases. 

An intriguing feature of the ES in phases with a broken con- 
tinuous symmetry is the relation [6] between its lower part 
and the so-called tower of states (TOS) energy spectrum ob- 
tained when a system with continuous symmetry breaking in 
the thermodynamic limit is placed in a finite volume [7]. We 
demonstrate that in the superfluid phase the ES envelope, sep- 
arated from the rest of the spectrum by an apparent "entan- 
glement gap", reflects the appropriate TOS Hamiltonian. This 
correspondence provides a framework to understand generic 
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FIG. 1: ES in the Mott phase, (left) Perturbative structure, DMRG data for equally bipartitioned cylinder, W = 8,L ^ 32, U ^ 100. 
Multiplicities are indicated near each group of states The dotted lines highlight the ES envelope, (center bottom) Bosonic configurations 
corresponding to leading order configurations for some ES levels. Only occupancies on sites nearer the boundary are shown (omitted sites 
are singly occupied). We restrict to SNa > 0. The arrows represent hopping (perturbative) processes to create a given configuration. In (/) 
we show the effective single particle (boundary) excitation corresponding to multiplet (6). (right) Entanglement dispersions. Data points are 
DMRG (U — 100); black lines are perturbative results, (i) Lowest ES level as a function of W. (ii) Entanglement dispersion for the lowest 
ES multiplet at SNa = 1. k denotes the momentum along the boundary. 



features of the ES envelope. 

The Bose-Hubbard model — We consider the two di- 
mensional Bose-Hubbard model on a cylinder of length L 
(with open boundary conditions) and circumference W (with 
periodic boundary conditions). The Hamiltonian is H = 
- i^lbj + ^-c-) + ^ii^i - 1) where are bosonic 
operators, rii — blbi, and U is the on-site repulsion. We re- 
strict to the case of unit filling. To calculate the ES we divide 
the cylinder in two parts A and B of size Va and Vb respec- 
tively, with W the boundary length (Fig. 1 center top). The 
number of bosons Na{Nb) in A{B) is a good quantum num- 
ber for the ES and can be used to label the ES levels. It is con- 
venient to introduce 6Na = Na — Va measuring the excess of 
bosons in A with respect to unit filling. Since the system has 
translational invariance along the boundary, the corresponding 
momentum fc is a good quantum number. 

The ES in the Mott insulator — In the U ^ oo limit, the 
ground state is a product state with one boson per site. The 
ES has only one level, = 0. As in Ref. [5], the ES at large 
U can be constructed through boundary-linked perturbation 
theory, treating the hopping term Hp = — (^I^j + ^-C-) 
as perturbation. 

Fig. 1 (left) plots the ES {£} for a system at C/ = 100 
(DMRG data). We denote with a the perturbative order giv- 
ing the leading contribution to the ES levels. As in gapped ID 
systems [5], the separation (^ log?/) between consecutive a 
reflects the perturbative nature of the ES. In the center bot- 
tom we show the configurations giving dominant contribution 
to selected groups of ES levels. The A part of these config- 
urations also gives the dominant configuration in the entan- 
glement eigenfunctions (eigenfunctions of pa)- At each 5Na 
and a > there is more than one level. The multiplicities are 



determined by the boundary length W, and can be understood 
in terms of hopping processes across the boundary. On the 
ES envelope (dotted line in Fig. 1 left), the group of levels at 
5Na is obtained by transferring 5Na bosons from subsystem 
B to A. This process appears at order 6Na (with amplitude 
jj-SNAy 'pjjg multiplicity ■m{SNA) is obtained as the number 
of ways of moving 6Na bosons across the boundary to give 
distinct configurations, i.e. m{6NA) — C{W,SNa) with C 
the binomial coefficient (see (a)(b) in Fig. 1 center). 

Similar boundary perturbative processes help explain ES 
levels above the envelope; c.f., (c)(d)(e) in Fig. 1 left and 
center. At a < 3 the leading order of all ES levels (not 
only the envelope) is given by pure boundary processes. At 
higher orders {a>3) processes further away from the bound- 
ary start contributing to the leading order of some ES levels, 
e.g., (e) in Fig. 1. The corresponding configurations involve 
hoppings perpendicular to the boundary, and thus are similar 
to the physics of ES configurations in ID chains [5]. This is 
reflected in the linear multiplicity (^W) in (e). Generically, 
the multiplicities grow exponentially with a {^W"). This 
growth of multiplicity is related to the area law of entangle- 
ment entropy which limits the performance of DMRG in 2D 
systems even in gapped phases. 

Entanglement dispersions & entanglement Hamiltonian — 
We now examine the substructures superposed on the per- 
turbative hierarchy of ES levels. As the simplest example, 
we consider the lowest ES states at SNa = 1, marked (b) in 
Fig. 1. These W states are composed from particle-hole exci- 
tations across the boundary at the W possible positions along 
the boundary. Fig 1 (f) shows how this local excess/depletion 
of bosons can be regarded as a boundary degree of freedom. 
Since the momentum k along the boundary is a good quan- 
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FIG. 2: (a) ES in the superfluid phase (DMRG data at (7 = 2 
for size L/4: — W). (b) The same ES as in (a) but after subtract- 
ing the contribution of the envelope, (c) Spacing 5 (defined in text 
and Fig. 3) plotted versus 1/W; various U. Dotted lines are fits to 
A/W + B/W^. (d) S as function of U. Dotted lines are fits to 
S ~ \/U (expected for 



turn number, in the entanglement eigenstates these excitations 
will not be localized, but will appear as momentum eigen- 
states, similar to a single-particle dispersion. This "entangle- 
ment dispersion" is given perturbatively up to 0{U~^) as 

= log(C/ V2) + 2{2W + 10)/[/2 - 24 [l /U + {2W- 

25)/[/^j cosfc-68/[/2cos2fc + 552/[/3cos3fc. (2) 

Comparison with DMRG data (U = 100, W = 8) is shown 
in Fig. 1 right panel (ii). 

The spectrum (2) can be interpreted as the spectrum of the 
"Hamiltonian" 'Hb = const. -\-^^^ Ar{a\ai+r + h.c.) where 
i labels sites of the boundary and a; are the boundary degrees 
of freedom. This entanglement Hamiltonian is a boundary 
Hamiltonian of single particle tight-binding form, with hop- 
ping amplitudes Ar decaying exponentially with distance r 
[8, 9]. These results can be generalized to each group of states 
on the ES envelope: the envelope states in sector 5Na are 
described by a boundary entanglement Hamiltonian T-Le for 
5Na particles. Expressions for T-Le can be obtained pertur- 
batively [9]. Up to lowest subleading order, is of nearest 
neighbor tight-binding form (on a chain). 

As the Mott-superfluid transition is approached, T-Le be- 
comes more and more long range [8, 9], similar to the cor- 
relation length which diverges upon approaching the phase 
transition. 

The ES in the superfluid phase — In the superfluid phase 
{U < 16.739 [10]) the ES looks dramatically different 
(Fig. 2). There is a clearer separation between a low-lying 
"envelope" and the rest of the ES, but the envelope now has 
quadratic dependence on SN^, and there is only a single en- 
velope level at each SNa- These features are due to the fact 
that the underlying superfluid state has spontaneous breaking 
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FIG. 3: (top) Restructuring of ES across the Mott-superfluid tran- 
sition (DMRG data for W = 8 = L/4, equal bipartitions). The 
arrows explain the definitions of the quantities A and 5. (bottom) A 
and 5 plotted versus U for several values of the boundary length W. 
The vertical line denotes the critical value Uc- 



of U{1) symmetry in the thermodynamic limit. We can ex- 
plain some of these features through a correspondence with 
the tower of states spectrum [7], which is the physical low- 
energy spectrum obtained when a system with spontaneously 
broken continuous symmetry is placed in a finite volume. 
Since the finite-size ES is plotted against quantum numbers 
whose conservation is spontaneously broken only in the ther- 
modynamic limit, it is naturally related to the TOS spectrum. 
The precise relation [6] is that for the lower part of the ES 
He ^ T~Lt/Te where "Ht is the TOS Hamiltonian, Te is an 
effective temperature given by Te = Vg/C with Vg the ve- 
locity of the gapless excitations (for the Bose Hubbard this 
is the sound velocity), and C ^ W is the linear size of the 
system. The form of Te reflects the finite size behavior of the 
sound-wave gap. For the Bose-Hubbard the TOS Hamiltonian 
is Ht ^ (<5^)^/(x^) where N is the total particle number 
operator and % = dn/d/z is the compressiblity. This suggests 
that He - ^{SNa? ~ iSNA)VivsXW). 

We characterize this scenario through the quantities A 
("TOS gap") and 6 ("envelope curvature"), defined pictorially 
in Fig. 3. Formally, S = ^k=o - S.o and A = ^k=2w/w - Cfc=o, 
with (,k=n, ^k=2Tilw in the lowest ES multiplet at a = 1, 
while ^0 is the lowest ES level. In the thermodynamic limit 
He - 0{1/W), then S - 0{1/W). Furthermore, in the 
limit [/ — )• at fixed W, using ~ U and Vg ~ VU 
one obtains S ^ VU. Finally, assuming that the lowest ex- 
citations above the ES envelope are sound-wave like one can 
write He ^ [Ht +Hsw]/Te where Hgw describes the sound 
wave excitations. Since Te ~ 1/W this suggests that the gap 
A remains finite in the thermodynamic limit, although loga- 
rithmic decay cannot be ruled out [11]. 

Numerical data shown in Figs. 2 and 3 supports this pic- 
ture. The presence of a finite gap A for all the values of W 
and the quadratic behavior of the ES envelope provide some 
confirmation of the tower of states picture. Figs. 2 (c) and 2 
(d) show good agreement with the predictions 5 ^ 1/W and 
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FIG. 4: (a,b) Renyi entropy 5^^* (DMRG data). In (b) the data of 
(a) is rescaled by boundary length W, and the continuous line is the 
perturbative result, (c) Contributions to Trp\ from the envelope ES 
levels (rhombi) and from the higher levels (circles). Vertical lines in 
each panel marks the transition point. 



In Fig. 2 (b) and 3 we notice also that the gap A hardly 
changes with W (although a ^ 1/ log W type of decay cannot 
be ruled out with these sizes.) Surprisingly, A is also constant 
as a function of SNa- We also note that, above the gap, the ES 
levels possess further band-like structures with a band of W 
levels of width comparable to A, slightly but distinctly sepa- 
rated from higher levels. Field theory arguments involving the 
dynamics of sound waves indicate similar structures [11]. 

At the [/ = point, the ES can be calculated exactly [9, 
12]. There is a single ES level for each 6Na, i.e., only the 
envelope survives and the rest of the ES is pushed to infinity. 
The wavefunction for U = is an exact Bose condensate; 
the highly symmetric situation is similar to the ferromagnetic 
case [13]. 

The Mott-superfluid transition — Across the phase tran- 
sition {Uc ~ 16.739 [10]), A and 5 show "dual" behaviors 
(Fig. 3 bottom). In the Mott insulating phase, A ^ l/W"^ 
(from (2) one has W'^/S. ^ 7r^/6), while in the superfluid A 
converges to a possibly nonzero value or vanishes only loga- 
rithmically with W [11]. In contrast, the Ty-dependence of 5 
is: (5 ~ const, in the Mott insulator, 5 — in the superfluid. 
The observed finite size scaling behavior of 6 suggests that 
this quantity might allow a scaling collapse, similar to the ID 
example of Ref. [14]. For the 2D transverse field Ising model 
the closing of a quantity similar to 5 at the quantum critical 
point has also been reported [15]. 

Entanglement entropies — In the Mott phase all the en- 
tropies can be constructed perturbatively in 1/U. Perturba- 
tive considerations for entanglement entropies in 2D gapped 
phases also appear in Refs. [16-19]. We first consider the 
lowest ES level which is the single-copy entanglement. 
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entropy of index 7. Up to 0{U^^) we find ^0 = cW with 
c = 4/f/2 +488/C/^, i.e. the expected area law behavior. Fig. 
l(i) compares with DMRG data ?aU — 100. The dependence 
on W in perturbation theory is tricky [9]. Remarkably, up to 
a = A the only dependence is the linear area law and higher 
powers in W disappear in the entropy. The same is true for 



the von Neumann entropy, 



64 , 2 
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and the Renyi entropies with 7 = 2, 3, 

S^^'^/W « 8/[/2 + 968/;/^ Sf/W « 6/[/2 + lAM/U^. 

In the superfluid (0 < C/ < Uc) the entropies are still ex- 
pected to exhibit an area law, although, due to the symmetry 
breaking, subleading logarithmic corrections arise [6, 19, 20], 

(n) 

i.e. S\ ^ C1W+C2 log W. Moreover, atU — there is pure 
logarithmic behavior Sa ^ l/21og7yA [9]. This implies that 
the area law coefficient vanishes as U ^ and the entangle- 
ment entropies are dominated by the logarithmic corrections 
up to larger systems sizes as U decreases. 

The transition between these different regimes is illustrated 
in Fig. 4 for s'x^ . In Fig. 4 (a) we show 5*^^' for 3 < < 8 
as a function of U. A sharp change is visible at the phase tran- 
sition. The data divided by W (Fig. 4 (b)) shows that in the 
Mott phase the entropy follows an area law and the behavior is 

fully reproduced by the perturbative result within the expected 

(2) 

regime (continuous line). In the superfluid, S\ shows strong 
finite size effects. In Fig. 4 (c), we show the contribution to 
Trp^ coming from the ES envelope (rhombi) and from the ES 
levels above (circles). The contribution of the non-envelope 
levels is much lower, especially in the superfluid (about four 

order of magnitude). This means that most of the contribution 

(2) 

to S'X' comes from the envelope, which is responsible of the 
subleading logarithmic correction [6]. 

These results have implications for the range of applicabil- 
ity of DMRG calculations. The logarithmic contribution is 
due to the replication of the ES in every SNa sector, which 
would be suppressed if 6Na was not a conserved quantity, 
as happens in DMRG calculations where the U{1) symmetry 
is enforced to be broken. (See Refs. [21, 22] for analogous 
calculations in magnetic systems.) Our results thus provide a 
basis for a more quantitative explanation as to why explic- 
itly symmetry-broken DMRG or MPS outperforms simula- 
tions exploiting the associated conserved quantum numbers in 
situations where continuous symmetry breaking takes place. 

Discussion — We have presented a thorough analysis of 
the ES of the 2D Bose-Hubbard model, in particular highlight- 
ing very different structures of the ES in the Mott and super- 
fluid phases. In the Mott phase, the ES can be understood per- 
turbatively starting from the U 00 product state limit. The 
ES is organized in a perturbative hierarchy and within each 
group there are dispersion structures reflecting the freedom of 
excitations to move along the boundary while retaining the 
same perturbation order In the superfluid regime, the ES re- 
flects the spontaneous symmetry breaking of the underlying 
state, and is related to the tower of states known from finite- 
size studies of physical spectra. We believe the major features 
addressed in the two phases are generic for a wide range of 
gapped and gapless 2D systems. 



5 



The results of the present work open up a number of open 
questions and research avenues. First, our analysis points to 
an advantage of DMRG calculations with enforced symme- 
try breaking; exploration of such techniques are yet to be at- 
tempted in detail. Second, it would be interesting to see the 
TOS structure appear for gapless phases where the symmetry 
breaking is more complex than the U{1) symmetry broken in 
our case. The ES envelope would then be expected to have a 
nontrivial degeneracy structure, which would provide a more 
stringent demonstration of the correspondence between the 
ES and the TOS spectrum. Finally, the robust emergence of 
an area law in perturbation theory (cancellation of W" terms 
with a > 1) is not fully understood at this point, and deserves 
to be addressed in a more general setting. 

We acknowledge stimulating discussions with T. Grover, 
M. MetHtski, and R.R.P. Singh. The DMRG simulations have 
been performed on machines of the platform "Scientific com- 
puting" at the University of Irmsbruck - supported by the 
BMWF. 
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Supplementary Materials 



OVERVIEW 



g = . Up to fourth order in 1//7 the dominant eigen- 
value (corresponding to the lowest ES level ^o) is given by 

;^{o,o} ^ 1-4W /U"^ -{AS^W -^W"^) /U^ +0{U-'^) (S.l) 



In these Supplements, we provide additional details and 
information on several aspects of the ES of the 2D Bose- 
Hubbard model: 

• We provide details on the perturbative expansions of the 
entanglement entropies in the Mott phase given in the 
paper In particular we show that the expected area law 
behavior seems to be built into the structure of pertur- 
bation theory. 

• We provide perturbative expressions in the Mott phase, 
up to the first few orders, for the entanglement Hamil- 
tonian describing the levels of the ES envelope. We 
demonstrate that up to the first few orders the entan- 
glement Hamiltonian has the form of a tight-binding 
ID Hamiltonian for hard-core bosons, whose degrees 
of freedom are localized at the boundary between the 
two subsystems (boundary locality). By examining the 
sectors 5Na = 1,2, we show that while in the Mott 
phase the entanglement Hamiltonian is short range it 
becomes long range as the Mott superfluid transition is 
approached. 

• We discuss the U — point. We show that the ES con- 
tains only one level in each sector SNa for which we 
provide analytic expressions. The behaviors of the en- 
tanglement entropies and 'gaps' (A and S) are discussed 
and compared with the interacting case U^O. 



ENTANGLEMENT ENTROPIES IN MOTT PHASE 

As discussed in the main text, the entanglement entropies 
can be calculated perturbatively at large U. Here we provide 
expressions up to fourth order in 1/U. This requires reduced 
density matrx eigenvalues Xi (or ES levels ^i) up to order a — 
2 as defined in Figure 1 (left panel) of the main text. 

We find that all the ES levels depend only on the bound- 
ary length W and not on the transverse size L, reflecting the 
boundary local nature of the ES. The perturbative expressions 
for ES levels contain terms nonlinear in W. Remarkably, at 
the orders we have calculated, these cancel when one con- 
structs the Renyi and von Neumann entropies, giving rise to 
the expected area law (cx W) behavior. 



The reduced density matrix eigenvalues 



We denote the eigenvalues of pA at perturbative order a 



{a.SNA} 



and the corre- 



and particle number sector SNa by A 
sponding ES levels as ^j^"'^-^-^^ — _ ^{"."^^a}^ -pj^^ index i 
labels the levels within each group at order a and sector SNa- 



As may be expected from the boundary locality of the ES, 
the eigenvalue depends only on the boundary length W, and 
not on the cylinder length L in the direction transverse to the 
boundary. It is noteworthy that (S.l) contains terms nonlinear 
in W. 

Q = 1 . There are W levels at a = 1 in the sector SNa = 
±1. Up to fourth order in 1/U we find 



ai,±i} _ 2 



148 8 48 
224 



224 

-^cos2fcj (S2) 



with kj = 2Tr j/W (j = 1, 2 . . . , W^). 

a ~ 2 . Finally, we need to know the expansion 

of A{2.±2} 

and A^^'°^. At a = 2 order there are W{W—l)/2 eigenvalues 
each in sectors SNa — ±2 and W{W — 1) in sector SNa — 
(Figure 1 of main article). All these eigenvalues are degener- 
ate up to fourth order: 



^{2,±2} ^ ^{2,0} 



4 



(S3) 



Entanglement entropies 

Using the perturbative results (S.l), (S.2), (S.3), we derive 
the corresponding expressions for the entropies. 

The single copy entanglement. The single copy entangle- 
ment 5*^°°'' is obtained from (S.l) as 



(oo) 



S' 



which has the expansion 



1okA{°-°> 



(S.4) 



-,(oo) 
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c72 



,488 
IF 

-W^ — 
[/6 



W 



[/6 

3 16 



— +0{U-') (S.5) 



Since (S.l) is correct only up to terms 0{U^^), then in (S.5) 
only the terms up to fourth order are meaningful. It is remark- 
able that up to fourth order (S.5) contains only linear terms in 
W (area law), even though in (S.l) a quadratic term ^ 
is present. This term cancels when expanding the logarithm 
in (S.4). The ^ U^^ terms apparently violate the area law but 
are not meaningful at this order and have to be discarded. At 
present we are unaware of any proof that higher powers of W 
must vanish at every perturbative order 
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Renyi entropies. The expansion for the other Renyi en- 



tropies S*^^ (7 > 2) is obtained using 



"^A — 



logTrp;^. 



(S.6) 



For 7 = 2 we have 
Tip\ = l-W- 



W 



968 



32 

W^—^+OiU-"). (S.7) 



C/2 C/4 

Note again the presence of quadratic terms which arise from 
the muhipHcity ( 

in (S.6) and expanding the logarithm we get 

:(2 



W^) of Xf'^^^ and Af Plugging 



T.. 8 ^^968 

W 1- W 



o{u- 



(S.8) 



Again, we find the area law behavior up to meaningful order, 
due to the cancellation of terms as observed above for the 
single copy entanglement. Similar observations hold for the 
7 = 3 Renyi entropy: 



Trp; 



1 - W- 



12 



W 



2928 ^2g^0(t;-5), (S.9) 



C/2 [/4 

Von Neumann entropy. Up to order four the von Neumann 
entropy is 



Sa = -Trp^ log PA 



w 



i5e{-l,l} i=l 
W(W-l)/2 

E E Af'^>logAp^^> 



5G{-2,2} 



1=1 
1^(1^-1) 



E 



Af^°>logAf'''> 



(S.IO) 



Using (S.l), (S.2), (S.3), and the identity 

w 

E cos(27r/Wrg) = WSrfi mod w , 

q=l 

with r an integer, we obtain 
1 



Sa = 4:W 



{/2 



1 1 2 



16 , 2 
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-5\ 



(s.ii) 



As observed for the Renyi entropies, although (S.l), (S.2), 
(S.3) contain non linear terms in W, the von Neumann entropy 
exhibits the area law behavior. 



ENTANGLEMENT HAMILTONIAN FOR THE ES 
ENVELOPE 

On the ES envelope in the Mott phase (Figure 1 of main 
article), at each nonzero 6Na there is a cluster of states which 
are described by an entanglement Hamiltonian and show cor- 
responding dispersions. Below we give details for these en- 
tanglement Hamiltonians. 
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FIG. SI: Entanglement Hamiltonian for the envelope ES multiplet 
in the sector 5Na ~ 1. We plot the absolute value of the hopping 
amplitudes \Ar\ versus the hopping range r, on log-linear scale. The 
empty symbols are DMRG data in the Mott insulating phase (for a 
system of size W = L/2 — 8). Perturbation theory results are 
shown with crosses. 



The SNa — 1 sector. The entanglement Hamiltonian is 
given by the single particle hopping Hamiltonian 



-Hi 



h.c.) 



(S.12) 



where i = 1, 2, . . . , labels the sites near the boundary be- 
tween the two subsystems and r = 0, 1, 2 . . . is the hopping 
range. Up to 0{U^'^) the hopping amplitudes are found to be 

Aq = log(f7V2) + 2(2T4^ + 10)/C/2 + 0(f7^'') (S.13) 



Ax = -12 



l/?7 + (2VF-25)/(2[/3) +0{U-^) 

A2 = -M/U^ + 0{U-^) 
A3 EE 226/[/3 + 0{U-^) 

In Fig. SI we compare the theoretical predictions (S.13) with 
DMRG data. We consider a system of dimensions W = 
L/2 = 8 and f7 = 20,50,100. In Fig. SI we show \Ar\ 
(for r — 0,1,2,3) obtained from (S.13). We also plot the 
hopping amplitudes \Ar\ obtained by fitting the DMRG data 
to the single particle tight binding dispersion 



6» = Ao + 2 E ^, 



27r . 

J cos —mj 



(S.14) 



The perturbative expressions S.13 describe the DMRG data 
extremely well for U = 50, 100. Not surprisingly, the agree- 
ment is less perfect for U = 20 (the transition to the superfluid 
regime is at [/ w 17). The hopping amplitudes \Ar \ are seen 
to decay exponentially in the Mott phase with the hopping 
range r, but longer range hopping becomes more and more 
relevant as the superfluid transition is approached. 

The SNa — 2 sector We provide the first two non zero 
orders (in 1/U) for the entanglement Hamiltonian describ- 
ing the levels of the ES envelope at SNa = 2. The block 
in the reduced density matrix corresponding to this sector can 
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FIG. S2: Entanglement dispersion in the envelope ES levels at 
SNa = 2. ES levels ^ are plotted against the momentum k along 
the boundary direction. Filled symbols are DMRG data atU — 100 
for W = L/2 = 8. Empty symbols denote the theoretical predic- 
tion (S.18). 



Since in our DMRG implementation we did not have access to 
the momentum, we assigned the momentum labels such that 
the data matches the theoretical prediction. (The procedure 
is not unique because ES levels with different momenta can 
be degenerate, but this is not important for present purposes.) 
Apart from small deviations due to higher order corrections in 
1/U, DMRG data are in agreement with (S.18). 

The general case SNa > 2. The result (S.18) can be ex- 
tended to all the other clusters of levels in the ES envelope. 
The block in the entanglement Hamiltonian describing the en- 
velope states at any SNa is given (up to the first subleading 
order) by the the tight-binding Hamiltonian for \SNa\ hard- 
core bosons: 

nEiSNA)^2\SNA\\og^+n' + ... (S.19) 



V2 



with H' given again by (S.17). 



be written as 



Pa{SNa = 2) 



4 



w 



E 

1=1 



12 



'-(4 a, 



h.c. 



1 

(S.15) 

where i labels the sites of subsystem A near the boundary, 
and a| are hard-core boson creation operators. The dots de- 
note higher order corrections. The corresponding block in the 
entanglement Hamiltonian "Hb = — log pA is 



neiSNA^ 2) =Uog 



U 

7i 



n' 



with 



12 



w 

E 

i=l 



(a|a 



h.c. 



(S.16) 



(S.17) 



Apart from the constant shift 4 log (J7/ a/2) the entanglement 
Hamiltonian is tight-binding Hamiltonian for two hard-core 
bosons, (S.17). Thus the ES levels in the envelope multiplet 
at SNa = 2 can be obtained from the energy spectrum of two 
free hard-core bosons on a periodic chain. The single particle 
spectrum of (S.17) is given by 



24 



(/+ i) if W even 



= 

with ; = 0, 1 . . 
is obtained as 



if W odd 



,W -1. The ES, i.e., the spectrum of (S.16) 

Alog{U/V2) + ei,+ei, 



(S.18) 



withZi <l2=0,l,---,W ~1. 

In Fig. S2 we show the ES levels in the envelope multiplet at 
SNa = 2 (DMRG data for U = 100 and Vt^ = L/2 = 8). We 
also show the theoretical result (S.18) plotting the ES versus 
the momentum along the boundary k = 2ttIi/W + 2tt12/W. 



ES AT THE U = POINT 

Given a system of N non interacting bosons (U — 0), the 
ground state is an exact Bose condensate. All the bosons are 
in the single particle state 



V 

E' 



Mm 



(S.20) 



with V — LW the total number of sites of the system. (Since 
we consider unit filling we also have TV = V.) Note that the 
"usual" Bose condensate in the state at zero momentum (fc = 
0) corresponds to the choice c,; = Vi in (S.20). This 

would give the correct ground state for non-interacting bosons 
on a torus. Since we are working on a cylindrical geometry 
this is not exactly correct in our case. 

The ES obtained from the Bose condensate in the 
state (S.20) contains only one level in each sector with fixed 
Na that is given by 



(Na) _ 



log 



:ei 



d^)"^(Ei-' 

(S.21) 

The ES (S.21) is continuously connected to the ES envelope 
upon switching on the interaction (at U ^ 0). 

In other words, in the limit [/ -> only the ES envelope 
survives and the higher parts of the ES are pushed to infinity. 
This impUes that the gap A diverges (A oo) in the U 
limit. In our DMRG data (Figure 3 bottom left in main 
article), we do not see a divergence for U as low as 2. This 
could be because the divergence starts at smaller U at these 
sizes, or that DMRG simulations with finite boson number 
cutoff per site is not able to capture this divergence. 

For the toric geometry (when q = l/VV ) the ES (S.21) 
becomes 



Nb 



N\ 



NaINbI 



NA^/NE 
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log 
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(S.22) 
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For our cylindrical geometry one obtains (S.22) only if block 
A is half of the system. Since this is the bipartition that we 
considered in the paper (and we expect that the physics re- 
mains unchanged for different bipartitions) we can use (S.22). 
Note that in the limit SNa = Na- N/2 < N the ES is ex- 
actly parabolic 

C(-^^=4log(^)4(^iV.r. (S.23) 

Strikingly, (S.23) implies 5 ^ ^/V, whereas at [/ > 
the presence of sound waves (and of the ES levels above the 
envelope) implies 5 ^ 1/W ^ l/VV (as discussed in the 



main text). 

As a final remark, we mention that the entanglement en- 
tropy for non-interacting bosons does not obey the area law 
and in the limit — > oo is given by 

Sa-^IoeNa. (S.24) 

This is consistent with the physics that the envelope in the 
gapless phase gives a logarithmic contribution to the entan- 
glement entropies, while the area law is due to sound waves. 
Since sound waves are missing at t7 = 0, the coefficient of 
the area law term also vanishes in this limit. 



